Microscopic picture of aging in SiC>2 
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We investigate the aging dynamics of amorphous SiC>2 via molecular dynamics simulations of a 
quench from a high temperature T\ to a lower temperature Tf. We obtain a microscopic picture 
of aging dynamics by analyzing single particle trajectories, identifying jump events when a particle 
escapes the cage formed by its neighbors, and by determining how these jumps depend on the waiting 
time i w , the time elapsed since the temperature quench to Tf. We find that the only ^-dependent 
microscopic quantity is the number of jumping particles per unit time, which decreases with age. 
Similar to previous studies for fragile glass formers, we show here for the strong glass former SiC>2 
that neither the distribution of jump lengths nor the distribution of times spent in the cage are 
tw-dependent. We conclude that the microscopic aging dynamics is surprisingly similar for fragile 
and strong glass formers. 

PACS numbers: 61.20.Lc, 61.20.Ja, 64.70.ph, 61.43.Fs 
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If a system is quenched from a high temperature 7] to 
a lower temperature Tf below the glass transition, crys- 
tallization is avoided and a glass is formed. The result- 
ing out of equilibrium (aging) dynamics has been hothy 
debated for the last decades and remains unclear 0, y|. 
Most previous studies on the aging dynamics investigated 
quantities which are averages over all particles in the 
system, such as mean squared displacement, incoherent 
intermediate scattering function, dynamic susceptibility, 
and energy . On the other hand much less is known 
about single particle dynamics during agin g. For col- 
loids, Cianci et al. investigated the structure 1(1 U| and 
Yunker et al. fl2| focused on irreversible rearrangements 
as function of waiting time t w . Warren and Rottler used 
computer simulations to investigate single particle hop- 
ping events for a binary Lennard- Jones mixture without 
shear as well as for polymers with and without shear [l3r - 
[l5| . To gain a more complete picture of the microscopic 
processes during aging, we study single particle hopping 
(jump) events for the very different glass former Si02- 
Whereas the systems of Warren and Rottler are fragile 
glass formers, SiC>2 belongs to the class of strong glass 
formers 

We determine the number of jumping particles per unit 
time, the jump length, and the time spent in a cage for a 
wide range of waiting times t w and for several choices of 
T and Tf . To study the aging dynamics of amorphous sil- 
ica we carried out molecular dynamics (MD) simulations 
using the BKS potential [16[ for the particle interactions. 
Starting from 20 independent fully equilibrated configu- 
rations at high temperatures T g {5000 K, 3760 K}, the 
system is quenched instantaneously to lower tempera- 
tures T f e {2500 K, 2750 K, 3000 K, 3250 K}. To keep the 
temperature at Tf constant and to disturb the dynamics 
minimally, the Nose-Hoover thermostat was applied only 
for the first 0.33 ns (NVT), and the simulation was con- 



tinued in the NVE ensemble for 33 ns during which Tf 
stayed constant. For more information on details of the 
simulation see 0. 

We focus on the microscopic dynamics at the lower 
temperature Tf by analyzing the single particle tra- 
jectories r n (t). During the production runs at Tf we 
stored average positions f n (ti) and fluctuations a n {ti) — 

for each particle n at times U = 
correspond to averages over 



vl{ tl )-(v n {ti)Y 



I x (0.00327 ns). Here (. . 
3200 MD steps and 2000 MD steps for the NVT and NVE 
simulation runs respectively. We then use the resulting 
f n (ii) to identify jump events. For example Fig. [1] shows 
the y-component of r n (ti) for n = 315; rectangular boxes 
indicate identified jumps. We define a particle n to un- 
dergo a jump if its change in average position 



satisfies 



Ar„ = \r n {ti) - r„(ti_ 4 )| 



Ar n > 3cr Q 



(1) 



(2) 



where a a is the average fluctuation size for particle type 
a G {Si,0}. Since a a is intended to be a measure of aver- 
age fluctuations during each particles rattling within its 
cage of neighbors, we first determine the estimate <jg St a 
by averaging (a n (ti)) 2 over all times ti of a given simu- 
lation run at Tf and over all particles of the same type 
a. We then determine a a by redoing the average over 
(cn(ij)) , but by averaging only over times for which 
(c n (ti)) < 3cr 2 st a which roughly excludes jumps from 
the average. Note that the definition of Eq . (GO ) is simi- 
lar, but not identical to our analysis in [13, EJ To ver- 
ify that our results are independent of the details of the 
jump definition, we replaced Eq. ([2]) with Ar n > V2a a 
and found indeed qualitatively the same results as they 
are presented here, for which we used Eq. ((2). 
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FIG. 1: (Color online) As an example for the time-averaged 
trajectory r^(ti) we show here the z-component r^7 for the 
oxygen atom n = 315 for a single simulation run at Tf = 2500 K 
which had been quenched from T = 3760 K. For clarity, only a 
fraction of the simulation time is shown. 
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FIG. 2: (Color online) Number of jumping particles N p per time 
At w as function of waiting time t w for the case of O-atoms 
and Ti = 5000 K (bold lines and symbols) and T = 3760 K 
(dashed thin lines). To be able to include on the logarithmic 
scale the data-point for the first time window at t w = 0, we 
plot -£f-(ty> = 0) instead at t w = 0.005 ns. For comparison the 
arrows indicate the equilibrium times (te in @]). 



We thus identify for all simulation runs all jump events 
occurring during the production run at Tf . For each jump 
event k we determine the particle rik jumping from aver- 
age position (r^T) 1 at time t k to average position (r„ fe ) f 
at time t k (see in Fig. [TJdark green and cyan circles). 

Our focus is on the dynamics of the system as it is 
aging over time. We investigate it via the jump events 
and their dependence on the waiting time t w , i.e. the 
time elapsed since the temperature quench to Tf. We 
divide the simulation run into waiting time windows, as 
indicated in Fig. [T] [21j . For each jump event k with jump 
time t\ we determine the waiting time window which 
includes t k (in Fig.Q]the light green waiting time window) 
and assign to this waiting time window the waiting time 
t w of the left border of the selected time window (in Fig.Q] 
red arrow). 

We therefore obtain jump statistics for each waiting 
time window starting at time t w and of duration At w (see 
Fig.[T]). In Fig. [2] we show the number of distinct particles 
jumpingper observation time At w as function of waiting 
time f w [22j . We find for all investigated Tf and both T\ a 
clear ^-dependence. With increasing waiting time 
decreases following roughly a power law until equilibrium 
is reached and -&f-(tw) becomes independent of i w and 
T;. The power law exponents are approximately the same 
for O- and Si-atoms in the range [— 0.6/ns, — 0.3/ns]. As 
one might expect, the larger T f the more particles jump 
and the earlier the equilibrium time t 3 eq , i.e. the time 



when levels off. For comparison we include in Fig. [2] 
the equilibrium times t^ q determined via the intermediate 
incoherent scattering function C g (i w ,i w + i) (t 



c __ 

eq 



^23 in 



[8|). We find i J eq ~ ig q , i.e. agreement between the mi- 
croscopic equilibrium time i J eq (single particle jumps) and 
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FIG. 3: (Color online) Jump length (AR) (see Eq. (0) and 
Fig- [lj) as function of waiting time t w for the case of T = 5000 K 
and O-atoms (bold lines and symbols) and Si-atoms (dashed thin 
lines). Similar to Fig.[2]we plot (AR)(t w = 0) at t w = 0.005 ns. 



the macroscopic equilibrium (C q includes a particle 
average). 

Next we test whether the t w -dependence manifests it- 
self also in a microscopic length scale. As sketched in 
Fig. [U we define the jump length of event k of particle 



rife jumping at time tf k from (v nk Y to (r„ fc ) to be 



AR k = (f^ 



(3) 



Similar to above, we investigate the ^-dependence of 
(AR) by including in the average only events for which 
t k belong to the same waiting time window. The result- 
ing Fig. [3] shows that (AR) for oxygen atoms (solid thick 
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FIG. 4: (Color online) Distribution of the jump length P(AR) 
for the case of 7] = 5000 K, T f = 2500 K and for O-atoms and 
in the inset for Si-atoms. Different colors indicate waiting time 

tw ■ 



lines with symbols) is independent of i w (with the only 
exception of the first time- window) , and for silicon atoms 
(dashed thin lines) (AP) is only slightly independent. 
This is in stark contrast to -^f- of Fig. [5J which shows 
strong ^-dependence. The ^-independence of AR holds 
true even for the distribution P(AR), both for O- and for 
Si-atoms, as shown in Fig. @] for the case of 7] = 5000 K, 
Tf = 2500 K. We find similar results for all other in- 
vestigated T; and Tf. Consistent with Fig. [3j we find 
only independence for t w ^ 0.02 ns (which corresponds 
in an experiment to the undetectable instant of an in- 
finitely fast quench). For i w > 0.02 an additional peak 
occurs at AR ss which is mostly due to reversible jumps 
(as defined in 17]). Furthermore we find exponential 
tails P(AR) ~ exp (— AP/Pd CC ay) with Pdccay ~ 0.8 and 
0.3 A for O- and Si-atoms respectively (similar to the re- 
sults for a binary Lennard Jones mixture [HI). 

With the conclusion from Figs. [3] and @] that the length 
scale AR is t w -independent, we investigate next the time 
scales associated with the single particle jumps. We de- 
fine the duration of a jump event k to be 



a*5 =A-A 



(4) 



(see Fig. [T]) and the time between successive jumps of the 
same particle 



— l k+l l k 



(5) 



that means the time spent in the cage before the same 
particle jumps again (see Fig. [1]). The resulting (Aid) 
and (Atb) are shown in Fig. [5j The time between jumps 
(Aib) is several magnitudes larger than (Aid). For com- 
parison with (Aib) we include arrows on the right to 
indicate ip q (i w = 23.98 ns) of 0, which is defined to 
be the time for which C 9 (i w ,i w -I- i^ q ) = 0.625. Since 



FIG. 5: (Color online) We show here average jump duration 
(Aid) (lower four curves) and time between successive jumps of 
the same particle (Aib) (top four curves) using the definitions of 
Eq. © and Eq. ([5|) and Fig. [TJ The arrows on the right indicate 
i? q (iw = 23.98 ns) of [8j|. We include Ai d (0 ns) and Ai b (0 ns) 
at t„ = 0.005 ns. 



(Aib) > t? q , we conclude that (Aib) is characterizing a- 
relaxation. As above, we determined the independence 
by averaging Ai *j and Ai^ for all jump events k for which 
t\ belongs to the same waiting time window. By choos- 
ing this definition of (Aib) we prevent artifacts due to 
the different time window sizes, because only t\ (instead 
of Ai^) is required to be in the time window of considera- 
tion. For large i w , however, the finite simulation run time 
itot = 33.33 ns, causes (Aib) to decrease for waiting times 
i w ~ (itot — Aib). Ignoring this itot-specific decrease, 
we therefore obtain the surprising result that (Aib) is 
independent of i w . This independence of i w holds not 
only for the average (Aib), but even for the whole dis- 
tribution P(Aib), as shown in Fig. [BJ Also in Fig. [5] we 
notice that P(At b ) ~ Ai" 1 at T t = 2500 K, whereas 
P(Ai b ) ~ exp (-Aib/idccay) at T f = 3250 K. In Fig. 
we show how P(Aib) plotted versus Aib changes with 
the final temperature, for a fixed i w = 8.75 ns. We ob- 
serve that at intermediate temperatures, i.e. Tf — 2750 K 
and Tf = 3000 K, there is a crossover from power law to 
exponential decay. For comparison we include in Fig. [7] 
the same arrows as in Fig. [21 which indicate the equilib- 
rium times i^j. The crossover time occurs approximately 

at the same time when xr"(*w) and C g (i w ,i w +t) reach 
equilibrium. A similar crossover has been observed for 
kinetically constrained models (see Fig. 10 of |19|) a nd 
for a binary Lennard- Jones mixture (see Fig. 2 of [20(). 

In summary, we obtain the following microscopic pic- 
ture of aging: both the distribution of jump length and 
the distribution of times spent in the cage P(Aib) are 
independent of waiting time i w (similar to the results 
of Warren and Rottler 13l4l5l|). Instead the only in- 
dependent microscopic quantity is the number of jump- 
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FIG. 6: (Color online) Distribution of times between jumps 
P(Atb) for O-atoms, Ti = 5000 K and for T f = 2500 K and 
in the inset for Tf — 3250 K. Different symbols (and colors) 
correspond to different waiting times t w . 




FIG. 7: (Color online) P(At b ) for fixed i w = 8.75 ns, T; = 
5000 K and for O-atoms as log-log plot in the main figure 
and as log-lin plot in the inset. Different symbols (and colors) 
correspond to different final temperature Tf. Dashed lines are 
power law fits with exponents —1.0, —0.9, —0.6, —0.3 and dot- 
dashed lines are exponential fits P(Atb) ~ exp (— Afb/tdecay) 
with f dccay = 10, 6, 2, 0.5 ns for T t = 2500, 2750, 3000, 3250 K 
respectively. As in Fig. [2] we include for comparison ar- 
rows which indicate the equilibrium times £^ [f§]. For clarity, 
P(Atb) has been shifted by a factor of 10 _1 /10" 3 /10" 6 for 
T f = 2750/3000/3250 K respectively. 



ing particles per time, which decreases with increasing 

[Hj]). This is 



t w (similar to the results of Yunker et al 
consistent with the first hop time results reported in [13- 
HEj ]. In agreement with kinetically constrained models 
P(A^b) shows a crossover from power law to exponential 
decay [19j. Our results for the strong glass former SiC>2 
are surprisingly similar to the fragile glass former results 
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